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Abstract 

Parity violating (PV) effects in neutron-deuteron radiative capture are studied using Desplan- 
ques, Donoghue, and Holstein (DDH) and effective field theory weak potentials. The values of PV 
effects are calculated using wave functions, obtained by solving three-body Faddeev equations in 
configuration space for phenomenological strong potentials. The relations between physical ob- 
servables and low-energy constants are presented, and dependencies of the calculated PV effects 
on strong and weak potentials are discussed. The presented analysis shows the possible reason for 
the existing discrepancy in PV nuclear data analysis using the DDH approach and reveals a new 
opportunity to study short range interactions in nuclei. 
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I. INTRODUCTION 



Low energy parity violating (PV) effects play an important role in understanding the main 
features of the Standard model. Many nuclear PV effects were measured and calculated 
during the last several years. Despite the fact that existing calculations of nuclear PV 
effects are in a reasonably good agreement with the measured ones, lately it became clear 
(see, for example [l-4] and references therein) that it is rather difficult to describe the 
available experimental data with the same set of weak nucleon coupling constants using the 
traditional DDH [5] weak meson exchange potential. 

As a possible solution for this problem, a new approach, based on the effective field theory 



(EFT), has 



papers l|, |4 



aeen introduced to parameterize the PV effects in a model independent way (see, 
Is] and references therein). The main goal of the EFT approach is to describe 
a large number of PV effects in terms of a small number of constants (LEC), which are the 
free parameters of the theory. Unfortunately, since the number of experimentally measured 
(and independent in terms of unknown LECs) PV effects in two body systems is not enough 
to constrain all LECs jz-lO]. In order to determine these constants it is necessary to include 
also the data obtained on heavier nuclear systems. 

Furthermore one should better understand PV effects in heavier nuclei because these 
effects might be essentially enhanced llMl3| in many body systems. However, how to apply 



the EFT approach for the calculations of PV effects in nuclei it is still an open question. 

To verify the possible issues related to the application of the DDH description of PV effects 
in nuclei and the possibility of systematic calculations of PV effects in nuclei using EFT 
approach, it is desirable to start from the calculations of PV effects in the simplest nuclear 
systems, such as neutron-deuteron (n-d) compound. PV effects for the elastic n-d scattering 
have been calculated recently 

HQ 

using both DDH and EFT approaches. However, 
before extending these techniques to the many-body nuclear systems, it is important to 
consider inelastic processes which are usually more sensitive to short range interactions. 



With this aim, we present in 



deuteron radiative capture 



16 



;his paper a comprehensive analysis of PV effects in neutron- 



191 ] using weak potential of DDH-type, as well as weak po- 



tentials obtained in pionless and pionful EFT. For the strong interaction, we have tested 
several realistic nucleon-nucleon potentials, also in conjunction with three-nucleon force. 
Three-nucleon wave functions have been obtained by solving Faddeev equations in configu- 
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ration space for the complete Hamiltonians comprising both weak and strong interactions. 

The paper is structured as follows. In the next section, a brief description of the employed 
formalism is presented. Then, we discuss the results of our calculations and perform a 
detailed analysis of model and cutoff dependence of the calculated PV parameters. In 
conclusion, the implications of our result are summarized. 



II. FORMALISM 



We consider three parity violating observables in the radiative neutron capture on 
deuterons (n + d — > 3 H + 7): circular polarization of the emitted photons (-P 7 ), asymmetry 
of the photons in relation to neutron polarization (a 7 ), and asymmetry of the photons in 
relation to deuteron polarization (AJ). For low energy neutrons, the expressions for these 
PV effects could be written in terms of parity conserving magnetic dipole (Ml) and parity 
violating electric dipole (El) transition matrix elements as: 
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Here the Ml and El amplitudes are defined as reduced matrix elements of the multipole 
operators 



XI j = {-q, J B \\fx\\J), with X = (M,E), 



(2) 



where and J are total angular momenta of bound state and scattering state respectively, 
and q is a momentum of the outgoing photon. The electromagnetic multipole operators in 
the limit of small q can be written as 



pMag 
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where J c (x) is a convection current, £i(x) is a magnetization current, p(x) is a charge op- 
erator, and q = u is the energy of photon. In our calculations, we use Ml operator up to 
N 3 LO in chiral order counting, which includes contributions from two-pion exchange and 



contact currents obtained in heavy baryon chiral perturbation theory [20j. For calculations 
of El amplitudes at the leading order, we use only El charge operator, which is related to 
3-vector currents by Siegert's theorem. Since, in the used spherical harmonics convention 
both parity conserving Ml and parity violating El amplitudes are purely imaginary, it is 
convenient to define real- valued Aij and Sj matrix elements as 

Ml j = *^=^=Mj, El j = -i^Sj, (3) 

VD7T\/47r V07T 



where fi N = 

The calculations of parity conserving Ml amplitudes for radiative n-d capture have been 



reported in papers 2(3, [2l| using the hybrid method, where the wave functions were obtained 
from phenomenological potential models and the current operators were derived from the 
heavy bary on chiral effective field theory. The results of these calculations can be approxi- 



mated 20j by the following expressions 1 



Mi = +21.87+ 10.76[(B model /B exp )- 2 - 5 - 1} fm*, 

Ma = -12.24 - 11.35[(B model /B exp )- 2 - 5 - 1} fm?, (4) 



where two low energy constants of two-body Ml operators are fixed [20| by experimental val- 
ues of 3 if and 3 He magnetic moments. In these expressions, Ml amplitudes and the binding 
energy of 3 H, B model depends on the strong interaction model. However, the observed ex- 
plicit correlation between the calculated The Ml amplitudes and the binding energy B mode i 
provides the unique opportunity to eliminate the model dependence. This might be done by 
setting in eq.fll]) B mode[ / B exp = 1. Ml amplitudes obtained in such a way lead to the value 
of the total neutron-deutron radiative capture cross section <j tot = 0.49(1) mb, which is well 
consistent with the experimental data. 

El amplitudes are calculated using three-body wave functions, which are obtained by 
solving Faddeev equations in configuration space. We have tested different combinations of 



1 The sign convention of Ml changed from [20( calculation to be consistent with the convention used in 
this work. 
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strong and weak potentials. For the strong (parity conserving) part of the Hamiltonian, we 
choose one of the realistic nucleon-nucleon interaction models, namely: AV18, Reid, Nijmll 
and INOY were employed. Also we have performed calculations for AV18 NN potential in 
conjunction with UIX three- nucleon force (denoted as AV18+UIX). For the parity violating 
part of the Hamiltonian - one of the weak potentials was employed, which was treated as 
perturbation. In this paper, we consider three types of parity violating weak potentials: the 
standard DDH potential with meson exchange nucleon-nucleon interactions, the potential 
derived from pionless version, and the potential derived from pionful version of effective 
field theory. Our approach could be considered as a hybrid method, which is similar to the 
hybrid approach in the line of Weinberg's scheme and which has been successfully applied for 



the calculations o 
hadronic systems 



weak and electromagnetic processes involving three-body and four-body 



20H25], as well as for calculations of parity violating 15[ and time reversal 



violating effects in elastic n-d scattering [26|, |27|]. It is worth mentioning that alternative 
calculations of parity violating effects in elastic n-d scattering using pionless EFT [28] are 
well consistent with the hybrid calculations [15], though the detailed comparison between 
these two methods is required. 



A. The parity violating potentials 



To understand the possible difference in the description of parity violating effects by 
DDH and EFT-type of potentials, we compare the operator structure of the potentials for 



the DDH potential 



29] and for two different choices of EFT potentials [1] which are derived 



from pionless and pionful EFT Lagrangian. All these potentials can be expanded in terms 



of Off operators 



14| as 



(n) 



a 



DDH, pionless EFT or pionful EFT 



(5) 



with the explicit forms for the operators Off and corresponding parameters c" listed in table 
H] 2 , where coefficients c" have dimension of [fm] and scalar functions have dimension 



Note that we changed the relation between the coefficient Cg in the weak Lagrangian [14| and the coef- 

^tCq, because of the inconsistency in 



151, c 



ficient c[ of weak potential from that of previous paper 

the convention. However, it does not affect our results in [15( because we calculated matrix elements of 
the operators O-, . 



5 



TABLE I: Parameters and operators of parity violating potentials, qa = 1.26, F w = 92.4 MeV. 
Tij = (3rf — Ti • Tj). Scalar function L\(r) = 3L\(r) — H\(r). 
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The operators in the last column are represented as products of isospin, spin, and 



vector operators X\j'±, which are defined as 



= [Pij, M r ij)] + i 



-(») 



«[Pij,/n(r<j)]-, 



(6) 



where Pij = ^zl. 

One can see that all weak potentials have the same structure, being represented by the 
15 basic operators which are allowed by symmetry. Thus, the difference between the weak 
potentials is due merely to the choice of coupling constants assigned to each operator and the 
scalar functions which describe the radial behavior of the term with the particular operator. 
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TABLE II: The cutoff parameters for the DDH parity violating potentials in GeV units 30]. For 
the masses of the mesons we use m n = 0.138 GeV, m p = 0.771 GeV, and = 0.783 GeV. 



DDH-I 
DDH-II 



K 

1.72 

oo 



1.31 

OO' 



1.50 

oo 



This leads to a principle difference between the DDH and EFT-type of weak potentials, 
because the scalar functions for the DDH potential are well constrained by a specific meson 
exchange forces, while the choice for the scalar functions in EFT approach is rather arbitrary. 
Therefore, EFT-potentials have more degrees of freedom, which leads to the possibility to 
parameterize a larger set of independent observables. However, one cannot predict the values 
of parity violating effects using EFT-potentials unless all LECs are determined. 

For the case of the DDH potential, radial functions f x {f), x = 7T,p, and u are usually 
written as normal Yukawa functions or modified Yukawa functions with corresponding cutoff 
terms 

^)^{e--e-[l + ^(l-f 

where, m x is a x-meson mass, and A^ is a corresponding cutoff parameter. We adopt two 
sets of the scalar functions, with and without cutoff terms, as described in Table HI1 and call 
them DDH-I and DDH-II. 

In the EFT, the results of the calculations of low energy observables should be indepen- 
dent of the specific form of the scalar functions /^(r) in the pionless EFT (/"EFT) potentials 
and of the form of the scalar functions used for the contact terms in pionful EFT (7rEFT), 
provided these functions are well localized (close to the delta function) and, at the same time, 
are smooth enough to be used in numerical calculations. This is because the dependencies 
on the mass scale (fi) and on the particular choice of the form of these functions must be 
absorbed by the renormalization of the low energy constants. Then, for our calculations in 
pionless EFT, we use two sets of the scalar functions, which we call //rEFT-I and /rEFT-II, 
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respectively 3 : 

Mr) = -^e->", for /EFT-I, 

A(r) = ^ A (r) = -J j^^e'^M AEFT-II (8) 

with mass scale parameters /x and A which provide a cutoff scale of the theory. For example, 
the natural scale of the cutoff parameters in pionless theory is (//, A) ~ m n . 

The pionful EFT model (7rEFT) has explicit long range interaction terms resulting from 
one pion exchange (VLi,lk) and higher order long range corrections (V^lr)- Also, it has 
middle range interactions due to the two pion exchange (V^mr), as well as a short range 
interactions {Vi^sr) due to nucleon contact terms. The radial part of the leading term of 
the long range one pion exchange, VLi^jj, is described by the modified Yukawa function 
f n (r). The short range interaction function V\ t sR m pionful theory has the same structure 
as pionless EFT. However, in spite of the structural similarity, the origins for these functions 
are different, therefore, as a consequence, their numerical values can be different. The only 
term in pionful EFT which has a different operator structure as compared to DDH or pionless 
EFT potentials, is a higher order long range correction term V^lr- However, we can ignore 
these higher order corrections from long range interactions, because they are suppressed and 
can be absorbed by renormalization of low energy constants [6|. Therefore, pionful EFT 
does not introduce a new operator structure as long as we neglect V^Yr term [f], The 
middle range interactions V^mr, or two pion exchange, can be described by functions L(q) 
and H(q) in momentum space 

\q\ \ 2m n J 4m^ + q z 

where, = (q°, q) = Pi — Pi = p£ — To transform these scalar functions into a 
configuration space representation by Fourier transform, we use the regulators S\(q). For 
the sake of simplicity, we use only one cutoff parameter for all regulators. Therefore, one 
can write 

1 r d 3 o A 2 
{L A (r), H A (r), f A (r), fl(r)} = —] j^e^ S A (q){L(q), H(q), 1, (10) 



Note that these terms are different from ones used in 



15] 



s 



where L\(r) and H\(r) correspond to two-pion exchange loop contributions, f\(r) and 
f][(r) describe short range contact terms and long range one-pion exchange contributions, 
correspondingly. It should be noted that we introduce the cutoff function even for the case of 
long range one-pion exchange potential to regularize a short range part of one-pion exchange. 
Among the possible choices, we use two types of regulators, which are called 7rEFT-I and 
vrEFT-II 4 : 

„2 (\2 _ Am 2 ) 2 

S^(q) = for ttEFT-I, Sf pole (q) = (A2 + ^ for ttEFT-II. (11) 

One can see that the function f\(r) in 7rEFT-I looks similar to the function for /7rEFT- 
II case; however, it leads to a different regularization since the typical value of the cutoff 
parameter for 7rEFT theory exceeds the pion mass scale and should be at least about of the 
p meson mass scale, while for pionless case it is close to the pion mass. Therefore, LECs for 
the same operators in pionless and pionful EFT potentials can be very different. 



B. Three nucleon wave functions 



Nuclear wave functions of initial (neutron-deuteron scattering) and final (bound triton) 
states of the neutron-deuteron radiative capture process are obtained in the context of 
non-relativistic quantum three particle problem. We consider neutrons and protons as the 
isospin degenerate states of the same particle nucleon, whose mass is fixed to h 2 /m = 
41.471 MeV-fm. The three-particle problem is formulated by means of Faddeev equations in 
configuration space . Using isospin formalism, three Faddeev equations become formally 
identical, which for pairwise interactions reads 

(E-H -V ij )i; k = V ij (^ l + ^), (12) 

where (ijk) are particle indexes, Hq is kinetic energy operator, V%j is a two body force 
between particles i, and j, ipk = 4>ij,k is so called Faddeev component. In the last equation, 
the potential formally contains both strong interaction, parity conserving, part (V? c ) and 
weak interaction, parity violating, part (V PV ), i.e.: Vij = V? c + V? v . Due to the presence 
of parity violating potential, the system's wave function does not have a definite parity and 



Note that the convention is different from the one used in 



15|- 
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contains both positive and negative parity components. As a consequence, the Faddeev 
components of the total wave function can be split in to the sum of positive and negative 
parity parts: 

i>k = ^t+^k ( 13 ) 

At low neutron energies, the dominant components of both initial and final state nuclear 
wave functions have positive parity. Parity violating interaction is weak {V^ v « V? c \ 
then by neglecting second order weak potential terms one obtains a system of two differential 
equations: 

(E-H - V? c ) = V? c (tf + (14) 
{E-Ho- V?°) r k = Vg°(i/>r + ^J) + V? v (tf + + ( 15 ) 



One can see that the first equation (fT4|) defines only the positive parity part of the wave 
function. This equation contains only a strong nuclear potential and corresponds to the stan- 
dard three nucleon problem: s-wave neutron-deuteron scattering, or an bound state of the 
triton. The solution of the second differential equation (fl5l) . which contains inhomogeneous 
term V? v (ipf + ijj^ + ip k ), gives us negative parity components of wave functions. 

To solve these equations numerically, we use our standard procedure, described in detail 
in [33j. Using a set of Jacobi coordinates, defined by x k = (rj —Vi) and y k = ^(tv 



we expand each Faddeev component of the wave function in bipolar harmonic basis: 

^ = E | 0* J,. H)« ® I fak tk )rr, • ( 16 ) 

where index a represents all allowed combinations of the quantum numbers presented in 
the brackets, l x and l y are the partial angular momenta associated with respective Jacobi 
coordinates, s« and U are spins and isospins of the individual particles. Functions F a (xk,yk) 
are called partial Faddeev amplitudes. It should be noted that the total angular momentum 
J, as well as its projection M, are conserved. Isospin breaking is taken fully into account 
by considering both T = 1/2 and T = 3/2 channels of the total isospin. 

Equations (fT4l) and (fi~5l) must be supplemented with the appropriate boundary conditions 
for Faddeev partial amplitudes F^. First of all, partial Faddeev amplitudes are regular at 
the origin: 

F±(0,y k ) = F±(x k ,0) = 0. (17) 
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For a bound state problem, system's wave function also vanish exponentially as either Xk or 
Uk becomes large. This condition is imposed by setting Faddeev amplitudes to vanish at the 
borders (x max ,y max ) of a chosen grid , i.e.: 

Fa ( x k, Umax) = 0, F±(x max ,y k ) = 0. (18) 

For neutron-deuteron scattering with energies below the break-up threshold, partial Faddeev 
amplitudes also vanish for — > oo, thus the last equality in ( |T8l) also applies for the 
scattering. 

At yk — > oo, all the Faddeev amplitudes vanish except for those consistent with the open 
channel, describing neutron-deuteron relative motion. For the case of thermal neutrons, we 
keep only the relative s-wave amplitudes in the asymptote. This behavior is imposed by: 

2 

F a ±} (x, y max ) = f!^l tSxttx (x)(y m ax - ■^=aj)Kfi 8 iv,^ 8 i^- ( 19 ) 



Here, fi^] x (%) are reduced deuteron wave function components with respective parity (±), 
orbital momenta l x , total angular momentum j x , total spin s x and total isospin t x . The 
corresponding deuteron wave function is calculated before three-nucleon scattering problem 
is undertaken. Neutron-deuteron scattering lengths aj for an angular momenta J = 1/2 
and J = 3/2 are obtained by solving equation (|T4l) . 

The formalism described above can be easily generalized to accommodate three-nucleon 



forces, as described in paper 



34]. 



C. Evaluation of the matrix elements 

In order to calculate the parity violating El matrix elements, we define real 8^ matrix 
elements corresponding to each operator as 

Sj = J2c n £ ( ;\ (20) 
/) 

where the sum is taken over different parity violating operators with corresponding LECs 
c n , defined in the table HI At the leading order, the electromagnetic charge operator does 
not violate parity. Therefore, the parity violating El amplitude results only from the small 
admixture of the parity violating component of wave functions. In the convention we use 
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the parity violating wave functions are purely imaginary both for bound state as well as for 
zero energy n-d scattering, one has 



£ 



(n) _ _p{n) , c(n) 
J ~ c J,(+)-|- c J,(-)> 



(21) 



where £/,(±) are amplitudes for transitions from a parity conserving scattering wave to a par- 
ity violating bound state, and from a parity violating scattering wave to a parity conserving 
bound state, respectively. 

In the first order of perturbation, parity violating El amplitudes can be presented as 
a linear combination of matrix elements calculated for each of the parity violating 

potential operators O^. Then, all PV observables a 7 , P 7 , Aj can be expanded as 

* = £ fe) *">. (*> 

(where X stands for a 7 , P 7 , or A^, and /ijv is introduced because of a dimension of coeffi- 
cients c m ) in terms of corresponding multipole amplitudes X^ m \ presented by the following 



expressions: 
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(25) 



(\M-^ + \Mf) 

It should be noted that for EFT potentials, each parity violating coefficient c n has an 
explicit cutoff or scale dependence multiplier -%( or -jg). Therefore, we present all results in 
normalized forms, as /i 2 (or A 2 ) x £( m )(or X^), to remove this artificial scale dependence. 

We calculate the parity violating El amplitude using 1-body charge operator 



(26) 



where, Qi and are i-th nucleons charge and position in the center of mass system, such 
that 

Q*»*< = 2 ( 2 a33 ^ 2 ~ T!^ 3 ^ 3 



(27) 



12 



Then, using the wave function expansion 



l^> = £^f^K (28) 



one obtains 



x||a)(/3|(r 2 -r 1 ) z |a) 



+ — <P||2/!I«><P| ( T3 „ — I a) , (29) 



a,/9 
1 

-=2, 
2V3 

where (^/f) 3 comes from the normalization of y. For these amplitudes, the integration 
over radial function is done numerically but angular parts of matrix elements are calculated 
analytically. 

III. RESULTS AND DISCUSSIONS 

The results of our calculations are presented separately for three choices of weak poten- 
tials: for the DDH potential, for the pionless and for the pionful potentials derived in the 
EFT approach. 

A. The DDH potential results 

The results obtained with the DDH potential are in a reasonably good agreement with the 



previous calculations [16l-ll9|. considering the difference in wave functions, and give us the 
opportunity to estimate the values of all PV effects in terms of PV meson-nucleon coupling 
constants h as 

a n = 0.42/4 - 0.17/iJ + 0.085/iJ + 0.008/iJ - 0.238/i° + 0.086/4 - O.OIO/^ 1 = 4.11 x 1(T 7 (30) 
P 7 = -1.05/4 + °- 19/i °, - 0.096/4 - °- 018/ 4 + 0.28/i° - 0.046/4 + 0.023/iJ = -7.31 x 10(31) 
A} = -1.51/4 + 0.17/i° - 0.083/4 - 0-024/4 + 0.024/i° + 0.013/4 + 0.032/iJ = -9.05 x l(ffi2) 

The coefficients in these expressions are obtained using strong AV18+UIX and weak 
DDH-II potentials, while the final values of PV observables are given for the "best" val- 
ues of DDH coupling constants. The contributions of different PV operators to transition 
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amplitudes £j,(p), where (P) indicates the parity of the scattering waves, are shown in ta- 
ble III II One can see that unlike the n-d elastic scattering case, there is no dominance of 
J = | channel and, as a consequence, all operators contribute almost equally to the capture 
process. 



TABLE III: Parity violating amplitudes £j^p) in fhi2 units, where (P) stands for the parity of the 
scattering wave, calculated with AV18+UIX strong and DDH-II weak potentials. 
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io- 2 


i 

i . 


ol 


X ru 


4 


1.03 


X 


10" 2 


-1.32 


X 


10" 2 


1.47 


X 


10~ 2 


-2. 


.87 


x 10~ 3 


5 


1.26 


X 


10~ 2 


-1.56 


X 


10~ 2 


1.75 


X 


10~ 2 


-3 


.79 


x 10~ 3 


6 


-2.03 


X 


10~ 3 


-8.85 


X 


10~ 3 


-1.85 


X 


IO" 3 


1. 


.51 


x IO" 3 


7 


-2.42 


X 


10~ 3 


-9.62 


X 


10~ 3 


-2.45 


X 


10~ 3 


1. 


.94 


x 10~ 3 


8 


-7.37 


X 


10~ 3 


2.43 


X 


10" 2 


-1.08 


X 


10~ 2 


9. 


.51 


x 10~ 3 


9 


-7.10 


X 


IO" 3 


1.24 


X 


io- 2 


-1.05 


X 


io- 2 


-2. 


.14 


x 10~ 3 


10 


9.79 


X 


IO" 3 


-1.25 


X 


io- 2 


1.39 


X 


io- 2 


-2. 


.71 


x IO" 3 


11 


1.20 


X 


10" 2 


-1.48 


X 


10~ 2 


1.67 


X 


10~ 2 


-3 


.61 


x 10~ 3 


12 


-2.75 


X 


10~ 3 


9.29 


X 


10~ 3 


-4.10 


X 


io~ 4 


-9. 


.10 


x 10~ 3 


13 


-3.05 


X 


10~ 3 


1.84 


X 


io- 2 


-1.96 


X 


10~ 3 


-1 


.53 


x IO" 2 
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TABLE IV: The DDH PV coupling constants in units of 10 7 [b! p contribution is neglected). Strong 

a 2 9 2 a 2 

interactions parameters are |^ = 13.9, = 0.84, |j£ = 20, k p = 3.7, and k u = 0. 

DDH Coupling DDH 'best' 4-parameter fit [25] 

hi +4.56 -0.456 

h° -11.4 -43.3 

h 2 p -9.5 37.1 

-1.9 13.7 
h x p -0.19 -0.19 

hi -1.14 -1.14 



TABLE V: Parity violating observables for different potential models with the DDH-best parameter 
values and Bowman's 4-parameter fits in 10 -7 units. 





DDH-best values 






4-parameter fits 






models 


a n 


P-y 


A d 


a n 


P 7 


A d 


AV18+ULX/DDH-I 


3.30 


-6.38 


-8.23 


1.97 


-2.16 


-1.81 


AV18/DDH-II 


4.61 


-8.30 


-10.3 


4.60 


-5.18 


-4.46 


AV18+UIX/DDH-II 


4.11 


-7.30 


-9.04 


4.14 


-4.71 


-4.09 


Reid/DDH-II 


4.74 


-8.45 


-10.4 


4.70 


-5.25 


-4.46 


NijmII/DDH-II 


4.71 


-8.45 


-10.5 


4.76 


-5.26 


-4.41 


INOY/DDH-II 


9.24 


-12.9 


-13.8 


17.5 


-17.9 


-13.5 



ts, we compare PV observables 



351 ] of weak coupling constants 



To check the possible model dependence of these resu 
for the "best" DDH values and for the 4-parameter fit 
(see table II VI) . For weak potentials, we used both DDH-I and DDH-II radial functions 
with strong interactions described by AV18, AV18+UIX, Reid, Nijmll, and INOY models. 
The results for these calculations are summarized in table |V] The difference in the values 
of PV effects for the "best" DDH values and for the 4-parameter fit proves that the PV 
effects in the radiative capture are very sensitive to the particular choice of the values of 
meson-nucleon coupling constants. We observed rather significant model dependence of the 
individual matrix elements. This model dependence indicates a possible serious problem in 
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the calculation of PV effects in nuclei, they require more thorough analysis of our calculations 
for EFT-type potentials presented in that follows. 

B. Pionless EFT potential results 

We start to analyze the EFT approach with the PV potentials obtained in pionless EFT by 
using scalar functions corresponding to two different schemes for cutoff procedure: /rEFT-I 
and /rEFT-II. Calculated PV amplitudes for these two weak EFT potentials for the same 
AV18+UIX strong interaction model are summarized in tables I VI I and IVHI The difference 
in /rEFT-I and ^rEFT-II results is not surprising because they have different forms of the 
scalar functions, the correct comparison of the results should be done for the products of 
these amplitudes with corresponding low energy constants. Then, the renormalization of the 
LECs can absorb the differences in the amplitudes. Unfortunately, we do not have enough 
experimental data to obtain these LECs. 



TABLE VI: Parity violating amplitudes £j^ P) for AV18+UIX stong interaction and PV /frEFT-I 
potential with [x = 138 MeV. 



op 










1 


-1.64 x KT 1 


-1.83 x KT 2 


-7.04 x 10~ 3 


-1.45 x 10" 1 


4 


2.68 x KT 1 


-2.74 x KT 1 


3.99 x lO" 1 


-9.64 x 10~ 2 


6 


-6.16 x 1(T 3 


-1.96 x KT 1 


-3.90 x 10~ 2 


5.30 x 10~ 2 


8 


-3.02 x KT 1 


4.07 x 10" 1 


-2.97 x 10" 1 


2.18 x 10" 1 


9 


-8.63 x 1(T 2 


1.74 x 1Q- 1 


-1.48 x 1Q- 1 


6.34 x 10~ 3 
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TABLE VII: Parity violating amplitudes £j^ P) for AV18+UIX stong interaction and PV /frEFT-II 
potential with fj, = 138 MeV. 



OD £ 1 i , \ 


<S 1 / \ 

i(-) 




f(-) 


1 -5.79 x KT 1 


-7.66 x KT 1 


4.09 x 10~ 2 


-1.36 x 10" 1 


4 4.28 x 10" 1 


4.62 x 10~ 2 


6.81 x 10" 1 


1.57 x 10" 1 


6 7.11 x 1(T 2 


-1.79 x KT 1 


-5.55 x 10~ 2 


-2.91 x 10~ 2 


8 -5.92 x KT 1 


6.57 x 10~ 3 


-5.15 x 10- 1 


6.43 x 10~ 2 


9 -1.45 x KT 1 


3.13 x KT 1 


-3.16 x 10" 1 


1.39 x lO" 1 


The contributions of different operators from these two weak EFT potentials with the 


same AV18+UIX strong potential to PV effects 


are shown in tables IVIIII and IIXI One can 


see that in the pionless EFT, all operators have approximately the same 


level of contribution 


to the PV effects, which is 


consistent with the results for the DDH mo 


del. 


TABLE VIII: Parity violatinj 


2; observables for AV18+UIX strong potential for /kEFT-I at \x = 138 


MeV. The results are in fm~ 


2 units. 






11 =T 


1 (n) 








2.17 x 1(T 2 


-5.52 x 10~ 2 


-7.93 x 10~ 2 


X 


-7.94 x 10~ 2 


6.55 x 10~ 2 


3.16 x 10~ 2 


X 


-2.81 x 1(T 2 


5.96 x 10~ 2 


8.01 x 10~ 2 


o 4mjv r<¥ 
J1 x 


1.04 x KT 1 


-1.03 x 10- 1 


-7.58 x 10~ 2 


X 


3.81 x 10~ 2 


-4.29 x 10~ 2 


-3.67 x 10- 2 
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TABLE IX: Parity violating observables for AV18+UIX strong potentials for /rEFT-II at A = 138 
MeV. The results are in fm~ 2 units. 



n ^ A 2 4 n) 

1 ~ 2 - 73 x 10_2 216 x 10_2 9 - 19 x 10_3 

4 ^.{Ct + Ct) -5.56 X 10- 2 2.19 xl0- 2 -2.32 x lO" 2 

6 -ins-Ct -3.69 x 10- 2 6.53 x 10" 2 8.08 x 10~ 2 

X 

8 -%^Cf 8.75 x 10~ 2 -6.76 x 10~ 2 -2.62 x 10~ 2 

X 

9 tffi-Cf 6.71 x 10~ 2 -5.02 x 10~ 2 -1.71 x 10~ 2 



C. Pionful EFT potential results 

The PV transition amplitudes calculated for strong AV18+UIX potential and PV pionful 
EFT potential with cutoff parameter A = 600MeV are presented in table |2D The results for 
PV observables are provided in tables IXII and Kill These tables reveal strong dependence 
on the choice of the scalar functions which, as was mentioned in the previous section, are 
expected to be absorbed by corresponding LECs. (For the comparison with pionless case, 
one shall take into account additional A 2 /A 2 multipliers in the coefficients of leading one-pion 
exchange operators which appear due to loop diagrams contributions in pionful EFT.) 
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TABLE X: El amplitudes calculated for AV18+UIX and vrEFT-I at A = 600 MeV in fma unit. 
(The A 2 multiplier is not included.) 



operator 










1 


-3.51 x 10" 1 


-7.40 x 10~ 2 


-1.15 x 10~ 2 


-2.86 x 10" 1 


4 


3.56 x 10^ 2 


-3.86 x 10~ 2 


4.97 x 10~ 2 


-7.90 x 10~ 3 


5 


3.43 x 10~ 2 


-4.36 x 10~ 2 


4.81 x 10~ 2 


-9.25 x 10~ 3 


6 


-7.37 x 10~ 3 


-3.18 x 10~ 2 


-6.42 x 10~ 3 


4.20 x 10~ 3 


8 


-2.65 x 10~ 2 


7.65 x 10~ 2 


-3.75 x 10~ 2 


2.84 x 10~ 2 


9 


-2.32 x 10~ 2 


4.48 x 10~ 2 


-3.51 x 10~ 2 


-5.06 x 10~ 3 


13 


-4.32 x 10~ 4 


6.26 x 10~ 2 


-6.46 x 10~ 3 


-4.36 x 10~ 2 


14 


-1.33 x 10~ 2 


5.33 x 10~ 2 


-6.11 x 10~ 3 


-4.66 x 10~ 2 


15 


1.27 x 10~ 2 


1.28 x 10- 1 


-2.19 x 10~ 2 


-9.03 x 10~ 2 



TABLE XI: PV observables for PV vrEFT-I potential and AV18+UIX strong potential at A = 600 
MeV. 



operator 


7(71) 

a n 


p( n ) 


d 


1 


4.12 x 10~ 2 


-1.06 x 10" 1 


-1.53 x lO" 1 


4 


-1.08 x 10~ 2 


1.03 x 10~ 2 


7.00 x 10~ 3 


5 


-1.14 x 10~ 2 


1.13 x 10~ 2 


8.12 x 10~ 3 


6 


-3.62 x 10~ 3 


7.51 x 10~ 3 


1.00 x 10~ 2 


8 


1.51 x 10~ 2 


-1.63 x 10~ 2 


-1.33 x 10~ 2 


9 


1.00 x 10~ 2 


-1.26 x 10~ 2 


-1.23 x 10~ 2 


13 


9.34 x 10~ 3 


-2.07 x 10~ 2 


-2.83 x 10~ 2 


14 


9.87 x 10~ 3 


-2.20 x 10~ 2 


-3.02 x 10~ 2 


15 


1.70 x 10~ 2 


-3.79 x 10~ 2 


-5.18 x 10~ 2 



19 



TABLE XII: PV observables for PV vrEFT-II potential and AV18+UIX strong potential at A = 600 
MeV. 



operator 


7(n) 

a n 


p(«) 
r i 




1 


2.10 x 10~ 2 


-5.62 x 10~ 2 


-8.20 x 10~ 2 


4 


-6.89 x 10~ 2 


6.53 x 10~ 2 


4.34 x 10~ 2 


5 


-6.44 x 10~ 2 


6.32 x 10~ 2 


4.46 x 10~ 2 


6 


-2.09 x 10~ 2 


4.47 x 10~ 2 


6.03 x 10~ 2 


8 


9.18 x 10~ 2 


-9.83 x 10~ 2 


-7.93 x 10~ 2 


9 


4.97 x 10~ 2 


-6.25 x 10~ 2 


-6.04 x 10~ 2 


13 


4.90 x 10~ 2 


-1.09 x 10- 1 


-1.49 x lO" 1 


14 


2.71 x 10~ 2 


-8.36 x 10~ 2 


-1.26 x 10" 1 


15 


1.10 x 10" 1 


-2.44 x 1Q- 1 


-3.33 x 10- 1 



D. Cutoff and model dependence 

The presented results reveal the model dependence of the calculated matrix elements, 
both on weak as well as on strong interaction. This model dependence has a different level 
of importance in calculating PV effects for different approaches. For the case of the DDH 
approach, the model dependence is directly related to the reliability of the calculations of 
PV effects in nuclei. In general, the EFT approach shall lead to model independent results; 
however, to guarantee the model independence, the intrinsic cutoff dependence must be 
checked explicitly. For the case of a "hybrid" EFT approach, which is not completely free 
from the possible model dependence, a careful analysis of both cutoff and model dependence 
of matrix elements and physical observables is required. 

In our approach, we used numerically exact wave functions of three-nucleon systems, 
however they depend on the choice of the strong Hamiltonian. Another possible source of 
the model dependence is the choice of PV violating potentials, which, for the EFT approach, 
means a choice of the scalar functions used for the regularization. It should be noted 
that in EFT, the model dependence of physical observables is not directly related to the 
model dependence of the calculated PV amplitudes because they are affected by the model 
dependence of the corresponding LECs. Unfortunately, at the present time these LECs are 
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unknown, which prevents derivation of PV observables. 

Since most phenomenological strong potentials have a similar long range behavior, cor- 
responding to one-pion exchange, the main difference between strong potentials is related 
to the middle and the short range contributions. Thus, rather strong model dependence 
of PV amplitudes implies that matrix elements related to n-d radiative capture process are 
sensitive to these short range interactions. This sensitivity to a short range dynamics is a 
new phenomenon observed in radiative n-d capture and is in direct contrast with the case of 
rarity violation in elastic n-d scattering where PV matrix elements are practically insensitive 



15| to the choice of the strong potential. 



This is partially related to the fact that in the case of the elastic n-d scattering, the 
dominant contribution to PV effects comes from the J = 3/2 channel, which is repulsive 
and thus less sensitive to the short range details of the potential. On contrary, in the case 
of n-d radiative capture, almost all channels contribute equally to the values of PV effects. 
In addition to that, for the radiative capture, the mechanism of pion exchange is not a 
dominant one, and, as a consequence, contributions from heavier meson exchanges (short 
distance contributions) become important. Therefore, one can see a number of reasons why 
PV three-body radiative capture processes should be more sensitive to the short distance 
dynamics than PV effects in three-body elastic scattering. It should be noted that even 
in the two-body case, a circular photon polarization P 7 in n-p radiative capture, which is 
not dominated by one-pion exchange, shows stronger model dependence 30] than a 7 , which 
relies on one-pion exchange contributions. 

As it is mentioned above, strong dependence of PV effects on the choice of potentials 
could be a serious problem in the case of the DDH meson exchange model, implying an 
uncertainty in the theoretical predictions and a difficulty in comparing results of different 
calculations. On the other hand, in a regular EFT approach, dependence on a cutoff and on 
the choice of a scalar function must be absorbed-copensated by the renormalization of the 
low energy constants. After the proper renormalization one must get model independent 
prediction of the low energy observables. This is not exactly true for the hybrid method, 
where strong interactions are introduced by a phenomenological strong potentials. However, 
it can be argued that the short distance details of the system dynamics would not be very 
important for the calculations of low energy observables according to the basic principle of 
the effective field theory. The removal of the model dependence, related to the difference in 
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short range parts of the wave functions, can be achieved by the introduction of the cutoff 
and renormalization of LECs in hybrid approach. A study of the behavior of the calculated 
matrix elements as a function of cutoff parameters in hybrid approach could be used to 
check the validity of these arguments. 




100 200 300 400 500 600 ' 100 200 300 400 500 600 

A(MeV) A(MeV) 

(a)/z 2 £ 3( + ) for operator 1 (b)^ 2 £ for operator 9 



FIG. 1: Cutoff and strong model dependencies of fj?£ in /m~2 for /rEFT-I calculated with AV18, 
AV18+UIX, Nijmegen-II, INOY, and Reid strong potentials. The multiplier fj? is used to absorb 
artificial cutoff dependence of c n coefficients. 

As an example, let us consider the fi 2 S3^ matrix elements as a function of a cutoff 
mass, which is calculated for operators 1 and 9 in the /ttEFT-I approach with different 
strong potentials (see Figfrj) . The choice of these operators is related to their symmetry 
properties: the operator 1 has quantum numbers corresponding to pion-exchange while the 
operator 9 - to p-meson exchange. It should be noted that since we use the same scalar 
functions both for the /rEFT-I and for the DDH-II schemes of calculations, we can apply 
the result of this analysis also to the calculations in the DDH-II scheme. Once again one 
observes rather strong dependence on the choice of a strong potential and on a cutoff mass 
parameter. 

Analyzing results of Fig JT] from the point of view of the DDH approach, where the matrix 
element for the operator 1 at p = corresponds to the pion-meson exchange and the matrix 
element for the operator 9 at p = m p corresponds to the rho-meson exchange, one can see 
a large strong potential model dependence for heavy meson exchange. This dependence 
indicates the importance of the inclusion of 3-body strong potentials. Unfortunately, most 
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calculations of PV effects in nuclear physics with the DDH potential do not include strong 3 



body forces, which could be a possible source for the existing discrepancy |36| in the analysis 
of PV effects. 

On the other hand, from the point of view of /rEFT, the reasonable cutoff mass scale 
cannot exceed the value of the pion mass, where the dependance on strong interaction 
potential is small. Since the cutoff in EFT could be considered as a measure of our knowledge 
of short range physics, increasing the cutoff parameter implies stronger dependence on the 
short distance details. Figfj] shows that by lowering the cutoff, one can diminish the strong 
potential model dependence. This is because by lowering of the cutoff parameter, we are 
effectively switching to the regime where the theory becomes sensitive only to a long range 
part of the interaction. Then, one can expect a smaller model dependence when the cutoff 
parameter is low, because all strong potentials have a similar long range behavior. Therefore, 
Figfj] is consistent with the basic principle of the EFT and shows that the hybrid method 
works well. 

The remaining weak dependence on strong interaction model at fi ~ m n scale could be 
related both to short and to long range parts of the potentials. If they are the remainder 
of the short distance part of wave function, the difference should be absorbed by LECs. On 
contrary the difference in the long range part of the wave function can not be removed by 
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the renormalization of LECs in the hybrid method. However, as demonstrated in 
long range part difference is governed by strong interaction observables and should be easily 
treated by analyzing the correlation between matrix elements and effective range parameters. 
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FIG. 2: Cutoff and strong model dependencies of A 2 £ in fm 2 for 7rEFT-I. The multiplier A 2 is 
used to absorb artificial cutoff dependence of c n coefficients. 



To analyze the possible model dependence for the pionful EFT approach, let us consider 
a contribution of operator 1 to A 2 £i ^ and to A 2 S 3 ^ calculated in 7rEFT-I approach 
(see Fig. |2]). In 7rEFT, the physical range for cutoff mass scale parameter A is about 
500 < A < 800 MeV. One can observe a rather important dependence on strong potential 
model in this region. We cannot discern long range from short range model dependency 
unless all LECs are determined. However, a smaller range of the variation of matrix elements 
for different strong potentials at the pion mass scale indicates that the contribution of the 
long range part of strong potentials to the region of the interest (500 < A < 800 MeV) is 
small. This means that the large model dependence in this range (500 < A < 800 MeV) is 
due to short range part of the wave function, therefore this cutoff and model dependence 
should be absorbed by higher order contact terms. 

Though the general behavior of the matrix elements are consistent with the expectations 
of EFT, the 3-body system is rather complicated one to see the direct relations between the 
2-body PV potential and 3-body PV matrix elements. Therefore, it is useful to re-analyze 
the two-body n-p capture process, for which the large model dependence for a circular 
polarization of photons, P 7 , was reported in 3o| . 
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E. Two body radiative capture(n + p — > d + 7) 



Parity violating asymmetry of photons for polarized neutron capture on protons and and 
their circular polarization for the case of unpolarized neutron capture can be written as 

-V2Re{M* 1 CS )E 1 ( 3 S 1 )} 

a ~ IM^W 2 1 } 

2Re[M 1 ( 1 S )P*( 1 S )] 

IM^So)! 2 ' 1 ] 

Here, we neglected M 1 ( 3 S'i), and i? 1 ( 1 P 1 ^— 3 Si) amplitudes. The Ei( 3 Si) amplitude is a 
sum of amplitudes with contributions from parity violating bound state wave function and 
from parity violating scattering wave ( 3 P\ ^— 3 Si). Since Mi( 3 Si) amplitude is suppressed, 
one can consider only i? 1 ( 1 S ) contribution to the P 7 , which is dominated by p and oj meson 
exchanges in the DDH formalism. (The a 7 is dominated by one-pion exchange.) 
The parity conserving Ml amplitude can be written as 

MlCS ) = i-^B—M = % ( ^^(393.06) fm'A . (35) 



Then, PV observables can be written as 



< - E fe) i-^r- 



E 

m 



5±) {-2yfo) E{ 3 1So) . (36) 



Using strong AV18 and weak DDH-II potentials, one can obtain 

a? n = 0.15^ + 0.00137^- 0.00405^-0.00137^, (37) 
= -0.0104/iJ- 0.00817^ + 0.0111/i°. (38) 
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TABLE XIII: Two-body Parity violating observables for potential models with DDH-best param- 
eter values and Bowman's 4-parameter fits. 



models oil P- 





DDH-best 


4-para. fit 


DDH-best 


4-para. fit 


AV18 +DDH-I 


5.25 x 1(T 8 


-4.91 x 10~ 9 


6.94 x 10~ 9 


4.76 x 1(T 9 


AV18 +DDH-II 


5.29 x 1(T 8 


-4.81 x 1(T 9 


1.76 x 1(T 8 


3.01 x 10~ 8 


NijmII+DDH-II 


5.37 x 1(T 8 


-4.99 x 1(T 9 


2.61 x 1(T 8 


6.41 x 10~ 8 


Reid+DDH-II 


5.33 x 1(T 8 


-4.85 x 1(T 9 


2.65 x 10~ 8 


4.68 x 10~ 8 


INOY+DDH-II 


5.60 x 1(T 8 


-3.94 x 1(T 9 


2.55 x 1(T 7 


9.68 x 10~ 7 



The calculated values of PV observables for different sets of strong potentials and different 
choices of DDH coupling constants are summarized in Table IXIIII One can see that the 
circular polarization P 7 , being dominated by heavy meson exchange, shows large model 
dependence in agreement with the analysis of n-d case. 
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FIG. 3: PV transition amplitudes /i 2 £j( + ) in fm 2 for AV18 strong potential and PV /kEFT-l 
potential. 

The cutoff and model dependence of the transition matrix elements calculated for oper- 
ators 1 and 9 using /rEFT-I approach, shown in Figure [31 remind the corresponding cutoff 
and model dependencies for the n-d capture process. One can see that the model depen- 
dence is more pronounced at the scale of p or u meson masses in comparison to the pion 
mass region scale. This is also consistent with the statement given for the hybrid approach 
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that one can regularize the short distance contributions by introducing a cutoff parameter 
and, as a consequence, reduce the uncertainty related to the short range interactions. This 
indicates that the possible reason for model and cutoff dependencies in 3-body (n-d capture) 
process has the same origin as those in 2-body case, and it can be treated regularly in the 
hybrid approach. 

For the completeness of the analysis, we present contributions of different PV operators 
to PV observables calculated in /rEFT and 7rEFT approaches with the AV18 potential (see 
tables IXIVI and [XV| correspondingly). The large difference between matrix elements with 
pionless and pionful PV potentials could be explained by different scales of cutoff parameters, 
and the comparison of th e results obtained from these two approaches could be done only 
after the renormalization of low energy constants. 



TABLE XIV: Two body Parity violating observables for AV18 and /frEFT potential at fj, = 138 
MeV. 



>rEFT-I 


yrEFT-n 


n 


2 (n) 
M On 




2 M 




1 


6.02 x 1(T 3 





1.36 x 1(T 2 





6 





2.48 x 1(T 2 





5.71 x 1(T 2 


8 





1.99 x 1(T 2 





1.78 x 1(T 2 


9 





-2.17 x 1(T 2 





-5.97 x 1(T 2 



TABLE XV: Two body Parity violating observables for AV18 and vrEFT potential at A = 600 
MeV. Only non-vanishing matrix elements are shown. 



op 


vrEFT-I 

7(n) 

a n 


ttEFT-II 

„t(") 


op 


ttEFT-I 

p(n) 


ttEFT-II 


1 


1.22 x 10~ 2 


7.14 x 10~ 3 


6 


3.05 x 10~ 3 


2.01 x 10" 3 


13 


1.11 x 10~ 3 


5.59 x 10~ 4 


8 


2.70 x 10~ 3 


1.80 x 10" 3 


14 


1.33 x 10~ 3 


7.41 x 10~ 4 


9 


-4.13 x 10~ 3 


-2.36 x 10~ 3 


15 


2.17 x 10~ 3 


1.00 x 10~ 3 
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IV. CONCLUSION 



PV effects in neutron-deuteron radiative capture are calculated for DDH-type and EFT- 
type, pionless and pionful, weak interaction potentials. Three-body problem was solved 
using Faddeev equations in configuration space, also by varying the strong interaction part 
of the Hamiltonian. Number of different phenomenological strong potentials has been tested, 
including AV18 NN interaction in conjunction with UIX 3-nucleon force. The analysis of the 
obtained results shows that the values of PV amplitudes depend both on the choice of the 
weak as well as strong interaction model. We demonstrated that this dependence has the 
expected behavior in the framework of the standard pionless and pionful EFT approaches. 
Therefore, this dependence is expected to be absorbed by LECs. Nevertheless, in order to 
obtain model independent EFT predictions for PV observables, one should perform all the 
calculations in a self-consistent way 37|. Using the "hybrid" approach we can minimize 
the model dependence, provided that all LECs are defined from the sufficiently large set of 
experimental data, which does not look practical in the nearest future. 

For the case of the DDH approach, the observed model dependence indicates intrinsic 
difficulty in the description of nuclear PV effects and could be the reason for the observed 
discrepancies in the nuclear PV data analysis (see, for example 38J and referencies therein). 
Thus, the DDH approach could be a reasonable approach for the parametrization and for 
the analysis of PV effects only if exactly the same strong and weak potentials are used in 
calculating all PV observables in all nuclei. However, the existing calculations of nuclear 
PV effects have been done using different potentials; therefore, strictly speaking, one cannot 
compare the existing results of these calculations among themselves. Further, most of the 
existing calculations do not include three body interaction which is shown to be important. 

We would like to mention that the observed sensitivity of PV effects to short range parts of 
interactions could be used as a new method for the study of short ranges nuclear forces. Once 
the theory PV effects is well understood, or once we use exactly the same parametrization 
for weak interactions, PV effects can be used to probe short distance dynamics of different 
nuclear systems described by different strong potentials. 
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